1 When does activity during the mating season start and end for each sex? (Figure 2, Top; Models 1-2)

1.1 Model 1: First, test effect of sex on date of first annual detection (i.e., start of swarm).

m1 = lmer(mindate2~sex + (1|site.wy2) + (1|pit_id),
              control=lmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)), 
              data = df.m1); summary(m1)
Linear mixed model fit by REML ['lmerMod']
Formula: mindate2 ~ sex + (1 | site.wy2) + (1 | pit_id)
   Data: df.m1
Control: lmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

REML criterion at convergence: 1119.1

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.5190 -0.5320 -0.0346  0.4860  3.9664 

Random effects:
 Groups   Name        Variance Std.Dev.
 pit_id   (Intercept) 0.19965  0.4468  
 site.wy2 (Intercept) 0.01754  0.1324  
 Residual             0.19144  0.4375  
Number of obs: 610, groups:  pit_id, 433; site.wy2, 7

Fixed effects:
            Estimate Std. Error t value
(Intercept)   1.0363     0.0599  17.301
sexF          0.2489     0.0766   3.249

Correlation of Fixed Effects:
     (Intr)
sexF -0.221
  • Result: Females start autumn activity at hibernacula later than males.

1.2 Model 2: Then, test effect of sex on date of last annual detection (i.e., end of swarm).

m2 = lmer(maxdate2~sex + (1|site.wy2) + (1|pit_id),
              control=lmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)),
              data = df.m2); summary(m2)
Linear mixed model fit by REML ['lmerMod']
Formula: maxdate2 ~ sex + (1 | site.wy2) + (1 | pit_id)
   Data: df.m2
Control: lmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

REML criterion at convergence: 1672.9

Scaled residuals: 
     Min       1Q   Median       3Q      Max 
-2.90593 -0.39386  0.05173  0.51876  2.39627 

Random effects:
 Groups   Name        Variance Std.Dev.
 pit_id   (Intercept) 0.15645  0.3955  
 site.wy2 (Intercept) 0.07805  0.2794  
 Residual             0.20263  0.4501  
Number of obs: 948, groups:  pit_id, 609; site.wy2, 10

Fixed effects:
            Estimate Std. Error t value
(Intercept)  2.71097    0.09313  29.111
sexF        -0.38722    0.06159  -6.287

Correlation of Fixed Effects:
     (Intr)
sexF -0.095
  • Result: Females end autumn activity at the hibernacula later than males.

1.3 Plot Figure 2, Top. Results of dates of mating activity analyses (Models 1-2)

  • Top panel plots data and model coefficients that determined the active period
  • Bottom panel visualizes peaks in activity by sex
print(p.fig2)

2 What factors contribute to sex differences in activity? (Figure 3A-C, Models 3-5)

2.1 Model 3, Figure 3A: First, what are the effects of temperature and sex on nightly activity?

# response: (1 = detected | 0 = undetected) 
m3 = glmer(detect2~bat.tavg*sex + (1|site.wy2) + (1|pit_id), family = binomial(),
            control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)),
            data = subset(df.m3)); summary(m3)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: detect2 ~ bat.tavg * sex + (1 | site.wy2) + (1 | pit_id)
   Data: subset(df.m3)
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
 22769.1  22816.7 -11378.5  22757.1    20773 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.1133 -0.5959 -0.4535  0.7844  4.1081 

Random effects:
 Groups   Name        Variance Std.Dev.
 pit_id   (Intercept) 0.6536   0.8085  
 site.wy2 (Intercept) 0.2089   0.4571  
Number of obs: 20779, groups:  pit_id, 544; site.wy2, 9

Fixed effects:
               Estimate Std. Error z value Pr(>|z|)    
(Intercept)   -0.687705   0.217297  -3.165  0.00155 ** 
bat.tavg       0.001265   0.007507   0.169  0.86618    
sexF          -3.759775   0.627498  -5.992 2.08e-09 ***
bat.tavg:sexF  0.150806   0.030333   4.972 6.64e-07 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) bt.tvg sexF  
bat.tavg    -0.662              
sexF        -0.169  0.233       
bat.tvg:sxF  0.160 -0.241 -0.980
# view temperature slopes for each sex using emtrends
kable(m3.emt)
sex bat.tavg.trend SE df asymp.LCL asymp.UCL
M 0.0012651 0.0075070 Inf -0.0134483 0.0159784
F 0.1520707 0.0294381 Inf 0.0943731 0.2097683
  • Result: Females are less active throughout their swarm period generally, but female activity increases with temperature. Male activity is less influenced by temperature.

2.2 Model 4, Figure 3B: Then, what are the effects of the length of the active period and sex on the probability of last detection?

# response: 1 = last day active | 0 = detected on subsequent nights
m4 = glmer(last.det2~isa*sex + (1|site.wy2) + (1|pit_id), family=binomial(),
           control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)),
           data=subset(df.m4));summary(m4)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: last.det2 ~ isa * sex + (1 | site.wy2) + (1 | pit_id)
   Data: subset(df.m4)
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2800.5   2839.4  -1394.3   2788.5     4787 

Scaled residuals: 
   Min     1Q Median     3Q    Max 
-1.665 -0.300 -0.142 -0.043 74.209 

Random effects:
 Groups   Name        Variance Std.Dev.
 pit_id   (Intercept) 5.469    2.339   
 site.wy2 (Intercept) 1.602    1.266   
Number of obs: 4793, groups:  pit_id, 433; site.wy2, 7

Fixed effects:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept) -4.29912    0.70400  -6.107 1.02e-09 ***
isa          0.70560    0.07336   9.618  < 2e-16 ***
sexM        -1.37247    0.52423  -2.618  0.00884 ** 
isa:sexM    -0.36702    0.06682  -5.492 3.97e-08 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
         (Intr) isa    sexM  
isa      -0.535              
sexM     -0.603  0.501       
isa:sexM  0.507 -0.950 -0.622
  • Result: Regardless of arrival date, females spend fewer days swarming than males before ending activity.

2.3 Model 5, Figure 3C: Last, what are the effects of temperature and sex on the probability of last detection?

# response: 1 = last day active | 0 = detected on subsequent nights
m5 = glmer(last.det2~bat.tavg*sex + (1|pit_id) + (1|site.wy2), family = binomial(),
            control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)),
            data = subset(df.m5)); summary(m5)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: last.det2 ~ bat.tavg * sex + (1 | pit_id) + (1 | site.wy2)
   Data: subset(df.m5)
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2647.9   2687.9  -1318.0   2635.9     5780 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-1.1904 -0.2993 -0.2391 -0.1375  9.1936 

Random effects:
 Groups   Name        Variance Std.Dev.
 pit_id   (Intercept) 0.01401  0.1184  
 site.wy2 (Intercept) 0.51982  0.7210  
Number of obs: 5786, groups:  pit_id, 542; site.wy2, 9

Fixed effects:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)    0.16202    0.56715   0.286  0.77512    
bat.tavg      -0.14453    0.02669  -5.416 6.11e-08 ***
sexF           4.67270    1.54263   3.029  0.00245 ** 
bat.tavg:sexF -0.16407    0.07701  -2.131  0.03313 *  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) bt.tvg sexF  
bat.tavg    -0.895              
sexF        -0.284  0.314       
bat.tvg:sxF  0.286 -0.321 -0.994
  • Result: Females will end activity on warmer nights compared to males that remain active even as temperatures decrease.

2.4 Plot Figure 3. Contribution of temperature and length of swarm on sex-specific activity (Models 3-5)

print(p.fig3)

3 How does mating phenology correspond to seasonal disease dynamics? (Figure 4A-B, Models 6-8)

3.1 Model 6: First, do pathogen loads differ between sexes during the autumn mating (active) season?

m6 = lmer(lgdL2~sex + (1|site.wy2), 
          data = df.m6); summary(m6) 
Linear mixed model fit by REML ['lmerMod']
Formula: lgdL2 ~ sex + (1 | site.wy2)
   Data: df.m6

REML criterion at convergence: 314.2

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-1.8127 -0.5967 -0.1576  0.5203  3.5268 

Random effects:
 Groups   Name        Variance Std.Dev.
 site.wy2 (Intercept) 0.2147   0.4634  
 Residual             0.3193   0.5650  
Number of obs: 172, groups:  site.wy2, 9

Fixed effects:
            Estimate Std. Error t value
(Intercept) -5.35210    0.17072 -31.350
sexF        -0.01397    0.09788  -0.143

Correlation of Fixed Effects:
     (Intr)
sexF -0.159
  • Result: No clear infection bias during autumn mating.

3.2 Model 7, Figure 4A: Then, how do pathogen loads change seasonally by sex?

m7 = lmer(lgdL2~sex*pdate2 + (1|site.wy2), 
          data = df2); summary(m7)
Linear mixed model fit by REML ['lmerMod']
Formula: lgdL2 ~ sex * pdate2 + (1 | site.wy2)
   Data: df2

REML criterion at convergence: 822

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.0626 -0.6530 -0.0952  0.6261  3.5533 

Random effects:
 Groups   Name        Variance Std.Dev.
 site.wy2 (Intercept) 0.2592   0.5091  
 Residual             0.5722   0.7564  
Number of obs: 343, groups:  site.wy2, 13

Fixed effects:
             Estimate Std. Error t value
(Intercept) -11.90983    0.69417 -17.157
sexM          3.06440    0.76088   4.027
pdate2        0.74914    0.06673  11.227
sexM:pdate2  -0.34035    0.07560  -4.502

Correlation of Fixed Effects:
            (Intr) sexM   pdate2
sexM        -0.754              
pdate2      -0.971  0.776       
sexM:pdate2  0.748 -0.992 -0.783
  • Result: Females develop more severe infections by early hibernation.

3.3 Model 8, Figure 4B: Last, does the end date of autumn activity influence infections in early hibernation?

m8 = lmer(lgdL2~med.maxdate + (1|site.wy2), 
           control=lmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)),
           data = df.m8); summary(m8)
Linear mixed model fit by REML ['lmerMod']
Formula: lgdL2 ~ med.maxdate + (1 | site.wy2)
   Data: df.m8
Control: lmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

REML criterion at convergence: 299.9

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.1009 -0.5816 -0.0953  0.6783  3.2574 

Random effects:
 Groups   Name        Variance Std.Dev.
 site.wy2 (Intercept) 0.6249   0.7905  
 Residual             0.6642   0.8150  
Number of obs: 114, groups:  site.wy2, 10

Fixed effects:
            Estimate Std. Error t value
(Intercept)  11.1252     3.9794   2.796
med.maxdate  -1.5740     0.4071  -3.867

Correlation of Fixed Effects:
            (Intr)
med.maxdate -0.998

3.4 Plot Figure 4. Seasonal sex-biased infection (Figure 4A-B; Models 7-8)

print(p.fig4)

4 Supplemental analyses and figures

4.1 Supp Figure 1: Visualize complete dataset

print(p.sf1)

  • Generally, males are more active than females throughout autumn mating, fully encompassing female swarm activity.

4.2 Use balanced dataset to support sex-specific effects

4.2.1 Supp Model 1, Supp Fig 2A: Test temperature-dependence of bat activity

sm1 = glmer(detect2~bat.tavg*sex + (1|site.wy2) + (1|pit_id), family = binomial(),
             control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)),
             data = subset(df1.tru)); summary(sm1)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: detect2 ~ bat.tavg * sex + (1 | site.wy2) + (1 | pit_id)
   Data: subset(df1.tru)
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4354.4   4392.6  -2171.2   4342.4     4266 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-1.5418 -0.5403 -0.4081 -0.2618  4.1241 

Random effects:
 Groups   Name        Variance Std.Dev.
 pit_id   (Intercept) 0.5480   0.7402  
 site.wy2 (Intercept) 0.0166   0.1288  
Number of obs: 4272, groups:  pit_id, 178; site.wy2, 8

Fixed effects:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)   -1.19533    0.43121  -2.772  0.00557 ** 
bat.tavg       0.02066    0.02126   0.972  0.33107    
sexF          -3.17607    0.74400  -4.269 1.96e-05 ***
bat.tavg:sexF  0.11847    0.03614   3.279  0.00104 ** 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) bt.tvg sexF  
bat.tavg    -0.964              
sexF        -0.573  0.565       
bat.tvg:sxF  0.563 -0.582 -0.980
#view temperature trends by sex
kable(sm1.emt)
sex bat.tavg.trend SE df asymp.LCL asymp.UCL
M 0.0206591 0.0212553 Inf -0.0210005 0.0623187
F 0.1391326 0.0293756 Inf 0.0815574 0.1967078

4.2.2 Supp Model 2, Supp Fig 2B: Differences in nightly temperature by sex and activity

sm2 = lmer(bat.tavg~detect2*sex + (1|pit_id) + (1|site.wy2), 
           data = subset(df1.tru)); summary(sm2)
Linear mixed model fit by REML ['lmerMod']
Formula: bat.tavg ~ detect2 * sex + (1 | pit_id) + (1 | site.wy2)
   Data: subset(df1.tru)

REML criterion at convergence: 19136

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.1646 -0.6908 -0.0348  0.7159  3.1546 

Random effects:
 Groups   Name        Variance Std.Dev.
 pit_id   (Intercept) 0.3333   0.5773  
 site.wy2 (Intercept) 1.3436   1.1591  
 Residual             4.9435   2.2234  
Number of obs: 4272, groups:  pit_id, 178; site.wy2, 8

Fixed effects:
              Estimate Std. Error t value
(Intercept)   19.35260    0.42333  45.715
detect21       0.13314    0.10616   1.254
sexF          -0.03922    0.12987  -0.302
detect21:sexF  0.48511    0.17724   2.737

Correlation of Fixed Effects:
            (Intr) dtct21 sexF  
detect21    -0.086              
sexF        -0.138  0.264       
detct21:sxF  0.051 -0.600 -0.322

4.2.3 Plot Supp Figure 2. Effects of temperature with balanced observations between sexes

print(p.sf2)

  • Results consistent with Model 3 in Figure 3A, using the balanced dataset between sexes:
    • Females increase nightly activity with temperature compared to males whose activity is relatively unaffected by temperature.
    • The difference between mean temperature on nights when bats were active compared to night when bats were undetected was greater for females than males. Females concentrate their activity to warmer nights compared to males that do not.

4.3 Does body condition affect nightly activity (Supp Figure 3; Supp Model 3)?

4.3.1 Test whether temperature and mass differently influence activity between sexes

sm3 <- glmmTMB(detect2~sex*mass*temp + (1|pit_id),
               family = binomial(),
               data = subset(df.sm3)); summary(sm3)
 Family: binomial  ( logit )
Formula:          detect2 ~ sex * mass * temp + (1 | pit_id)
Data: subset(df.sm3)

     AIC      BIC   logLik deviance df.resid 
  5146.4   5203.6  -2564.2   5128.4     4249 

Random effects:

Conditional model:
 Groups Name        Variance Std.Dev.
 pit_id (Intercept) 0.9313   0.965   
Number of obs: 4258, groups:  pit_id, 336

Conditional model:
                Estimate Std. Error z value Pr(>|z|)  
(Intercept)     -4.82732    2.28755  -2.110   0.0348 *
sexF           -19.80829    9.78416  -2.025   0.0429 *
mass             0.34114    0.25653   1.330   0.1836  
temp             0.10615    0.11711   0.906   0.3647  
sexF:mass        2.20187    1.17849   1.868   0.0617 .
sexF:temp        1.03982    0.49003   2.122   0.0338 *
mass:temp       -0.00289    0.01312  -0.220   0.8256  
sexF:mass:temp  -0.12246    0.05921  -2.068   0.0386 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

4.3.2 Plot Supp Figure 3: Effects of mass and temperature on nightly activity

print(p.sf3)

  • Results:
    • Females with lower body condition (ie body mass) concentrated activity to warm nights, compared to females with higher body condition.
    • The temperatures at which males were active was not affected by their body condition.
    • Females adjust nightly activity according to their body condition but males do not, suggesting the sexes are likely budgeting their energy differently.

4.4 Additional analysis showing association between between mating phenology and disease (Supp Figure 4; Supp Figure 4)

4.4.1 Test whether median end date of autumn activity at a site influenced mean pathogen loads in early hibernation

sm4 = lmer(mean.lgdL~med.maxdate + (1|site.wy2), data = subset(df2.sum)); summary(sm4)
Linear mixed model fit by REML ['lmerMod']
Formula: mean.lgdL ~ med.maxdate + (1 | site.wy2)
   Data: subset(df2.sum)

REML criterion at convergence: 31.1

Scaled residuals: 
     Min       1Q   Median       3Q      Max 
-0.82223 -0.53563 -0.06433  0.47551  1.18360 

Random effects:
 Groups   Name        Variance Std.Dev.
 site.wy2 (Intercept) 0.5788   0.7608  
 Residual             0.1053   0.3244  
Number of obs: 16, groups:  site.wy2, 10

Fixed effects:
            Estimate Std. Error t value
(Intercept)  11.0041     3.1930   3.446
med.maxdate  -1.5615     0.3285  -4.753

Correlation of Fixed Effects:
            (Intr)
med.maxdate -0.997

4.4.2 Plot Supp Figure 4. Effect of median end date of autumn activity on mean early hibernation infections at site.

print(p.sf4)

  • Results: Consistent with Figure 4B, infection severity in early hibernation is more severe at sites where autumn activity ended sooner.

4.5 How many nights of activity do individuals accumulate throughout the length of their active period?

4.5.1 Plot Supp Figure 5: Visualize cumulative number of active nights by sex

print(p.sf5)

4.5.2 Calculate average number of nights detected among individuals of each sex

kable(n.dates.sum)
sex mean_n.dates
F 3.810526
M 13.947471
  • Males continue to be active without any clear threshold in the number of dates active, whereas females cease activity after an average of 4 nights.

4.6 Model comparisons and validation

4.6.1 Compare activity period models as additive or interactive with site-year and sex as fixed effects (Supp Table 2)

First, begin date of activity models

m1b.add = lmer(mindate2~site.wy2 + sex + (1|pit_id), data = df.m1)
m1b.int = lmer(mindate2~site.wy2 * sex + (1|pit_id), data = df.m1)
  • Compare using AIC
kable(AIC(m1b.add,m1b.int))
df AIC
m1b.add 10 1142.960
m1b.int 16 1160.287

Then, end date of activity models

m2b.add = lmer(maxdate2~site.wy2 + sex + (1|pit_id), data = df.m2)
m2b.int = lmer(maxdate2~site.wy2 * sex + (1|pit_id), data = df.m2)
  • Compre using AIC
kable(AIC(m2b.add,m2b.int))
df AIC
m2b.add 13 1692.986
m2b.int 20 1711.309

4.6.2 View results of null model comparisons and AUC estimates from k-fold cross validation of binomial models

kable(nullcomps_full)
df AIC
m1.null 4 1134.2857
m1 5 1129.1469
m2.null 4 1715.0409
m2 5 1682.9043
m3.null 3 22821.9059
m3 6 22769.0835
m4.null 2 3623.1434
m4 6 2800.5278
m5.null 2 2741.8440
m5 6 2647.9294
m7.null 3 966.0277
m7 6 833.9688
m8.null 3 319.9266
m8 4 307.8518
  • All AICs of reported models are >2 scores below the null models (structured with response~1), indicating improvement over the null.

4.6.3 View AUC scores derived from k-fold cross validation of the logistic models reported in the main results.

kable(kfold_full)
Model AUC
m3 70.54%
m4 75.94%
m5 83.09%
  • Results indicate that >70%, 75%, and 83% of our test data was successfully predicted by the training models of Model 3, Model 4, and Model 5, respectively.

4.7 View sampling distributions by site, year and sex where applicable by model

ST3A. Model 1 Sample summary: Start dates of activity
Year Females Males
BA IT
2021 11 60
2022 16 134
MA N
2021 13 33
2022 10 72
NE MI
2021 2 36
2022 17 88
2023 26 92
ST3B. Model 2 Sample summary: End dates of activity
Year Females Males
BA IT
2020 4 63
2021 12 102
2022 17 139
MA N
2020
10
2021 18 73
2022 16 98
NE MI
2020
20
2021 7 73
2022 24 107
2023 35 130
ST3C. Model 3 Sample summary: Prob of nightly activity
Year Females Males
BA IT
2020 4 63
2021 12 101
2022 17 139
MA N
2020
10
2021 18 71
2022 16 100
NE MI
2020
20
2021 7 73
2022 24 107
ST3D. Model 4 Sample summary: Prob of last detection
Year Females Males
BA IT
2021 11 60
2022 16 134
MA N
2021 13 33
2022 10 73
NE MI
2021 2 36
2022 17 88
2023 26 92
ST3E. Model 5 Sample summary: Prob of last detection
Year Females Males
BA IT
2020 4 63
2021 12 97
2022 17 139
MA N
2020
10
2021 18 67
2022 16 100
NE MI
2020
20
2021 6 72
2022 24 107
ST3F. Model 7 Sample summary: Seasonal infections
Year Season Females Males
BA IT
2019 Early hiber 10 9
2020 Autumn mating 26 66
2020 Early hiber 5 12
2021 Autumn mating
3
2021 Early hiber 1 10
2022 Autumn mating 11 38
2022 Early hiber 3 13
MA N
2019 Early hiber 4 15
2020 Autumn mating 12 28
2020 Early hiber 5 11
2021 Autumn mating 14 12
2021 Early hiber 2 7
2022 Autumn mating 22 48
2022 Early hiber 4 9
NE MI
2019 Early hiber 5 13
2020 Autumn mating 5 37
2020 Early hiber 1 11
2021 Autumn mating
6
2021 Early hiber 5 8
2022 Autumn mating 23 34
2022 Early hiber 10 6
2023 Early hiber
20
ST3G. Model 8 Sample summary: Phenology-dependent early hibernation infections
Year N
BA IT
2020 16
2021 8
2022 10
MA N
2020 8
2021 9
2022 10
NE MI
2020 7
2021 12
2022 15
2023 19
LS0tCnRpdGxlOiAiQXBwZW5kaXggZm9yIE1hdGluZyBzeXN0ZW1zIGFuZCBzZXgtYmlhc2VkIGRpc2Vhc2UiCm91dHB1dDogCiAgaHRtbF9ub3RlYm9vazoKICAgICAgbnVtYmVyX3NlY3Rpb25zOiB5ZXMKICAgICAgdG9jOiB5ZXMKICAgICAgdG9jX2Zsb2F0OiB5ZXMKICAgICAgdGhlbWU6IHVuaXRlZAotLS0KCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFLCB3YXJuaW5nID0gRkFMU0V9Cm9wdHNfY2h1bmskc2V0KHdhcm5pbmcgPSBGQUxTRSwgbWVzc2FnZSA9IEZBTFNFKQpgYGAKCmBgYHs9aHRtbH0KPHN0eWxlIHR5cGU9InRleHQvY3NzIj4KICBib2R5ewogIGZvbnQtc2l6ZTogMTJwdDsKfQo8L3N0eWxlPgpgYGAKCiMgV2hlbiBkb2VzIGFjdGl2aXR5IGR1cmluZyB0aGUgbWF0aW5nIHNlYXNvbiBzdGFydCBhbmQgZW5kIGZvciBlYWNoIHNleD8gKEZpZ3VyZSAyLCBUb3A7IE1vZGVscyAxLTIpCiMjIE1vZGVsIDE6IEZpcnN0LCB0ZXN0IGVmZmVjdCBvZiBzZXggb24gKipfZGF0ZSBvZiBmaXJzdCBhbm51YWwgZGV0ZWN0aW9uXyoqIChpLmUuLCBzdGFydCBvZiBzd2FybSkuCgpgYGB7cn0KbTEgPSBsbWVyKG1pbmRhdGUyfnNleCArICgxfHNpdGUud3kyKSArICgxfHBpdF9pZCksCiAgICAgICAgICAgICAgY29udHJvbD1sbWVyQ29udHJvbChvcHRpbWl6ZXI9ImJvYnlxYSIsIG9wdEN0cmw9bGlzdChtYXhmdW49MTAwMDAwKSksIAogICAgICAgICAgICAgIGRhdGEgPSBkZi5tMSk7IHN1bW1hcnkobTEpCmBgYAoKLSAqKl9SZXN1bHRfKio6IF9GZW1hbGVzIHN0YXJ0IGF1dHVtbiBhY3Rpdml0eSBhdCBoaWJlcm5hY3VsYSBsYXRlciB0aGFuIG1hbGVzLl8KICAKIyMgTW9kZWwgMjogVGhlbiwgdGVzdCBlZmZlY3Qgb2Ygc2V4IG9uICoqX2RhdGUgb2YgbGFzdCBhbm51YWwgZGV0ZWN0aW9uXyoqIChpLmUuLCBlbmQgb2Ygc3dhcm0pLgoKYGBge3J9Cm0yID0gbG1lcihtYXhkYXRlMn5zZXggKyAoMXxzaXRlLnd5MikgKyAoMXxwaXRfaWQpLAogICAgICAgICAgICAgIGNvbnRyb2w9bG1lckNvbnRyb2wob3B0aW1pemVyPSJib2J5cWEiLCBvcHRDdHJsPWxpc3QobWF4ZnVuPTEwMDAwMCkpLAogICAgICAgICAgICAgIGRhdGEgPSBkZi5tMik7IHN1bW1hcnkobTIpCmBgYAoKLSAqKl9SZXN1bHRfKio6IF9GZW1hbGVzIGVuZCBhdXR1bW4gYWN0aXZpdHkgYXQgdGhlIGhpYmVybmFjdWxhIGxhdGVyIHRoYW4gbWFsZXMuXwoKIyMgUGxvdCBGaWd1cmUgMiwgVG9wLiBSZXN1bHRzIG9mIGRhdGVzIG9mIG1hdGluZyBhY3Rpdml0eSBhbmFseXNlcyAoTW9kZWxzIDEtMikKICAtIFRvcCBwYW5lbCBwbG90cyBkYXRhIGFuZCBtb2RlbCBjb2VmZmljaWVudHMgdGhhdCBkZXRlcm1pbmVkIHRoZSBhY3RpdmUgcGVyaW9kCiAgLSBCb3R0b20gcGFuZWwgdmlzdWFsaXplcyBwZWFrcyBpbiBhY3Rpdml0eSBieSBzZXgKYGBge3J9CnByaW50KHAuZmlnMikKYGBgCgo8ZGl2IGNsYXNzID0gImFsZXJ0IGFsZXJ0LWluZm8iIHJvbGU9ImFsZXJ0Ij4KICAtICoqX0tleSBmaW5kaW5nczpfKiogX0ZlbWFsZXMgaGF2ZSBhIHNob3J0ZXIgYWN0aXZlIHBlcmlvZCwgbWFsZSBhY3Rpdml0eSBmdWxseSBlbmNvbXBhc3NlcyBmZW1hbGUgYWN0aXZpdHksIGFuZCBtYWxlcyByZW1haW4gaGlnaGx5IGFjdGl2ZSBsYXRlciBpbnRvIGF1dHVtbi5fCjwvZGl2PgoKIyBXaGF0IGZhY3RvcnMgY29udHJpYnV0ZSB0byBzZXggZGlmZmVyZW5jZXMgaW4gYWN0aXZpdHk/IChGaWd1cmUgM0EtQywgTW9kZWxzIDMtNSkKCiMjIE1vZGVsIDMsIEZpZ3VyZSAzQTogRmlyc3QsIHdoYXQgYXJlIHRoZSBlZmZlY3RzIG9mIHRlbXBlcmF0dXJlIGFuZCBzZXggb24gKipfbmlnaHRseSBhY3Rpdml0eV8qKj8KCmBgYHtyfQojIHJlc3BvbnNlOiAoMSA9IGRldGVjdGVkIHwgMCA9IHVuZGV0ZWN0ZWQpIAptMyA9IGdsbWVyKGRldGVjdDJ+YmF0LnRhdmcqc2V4ICsgKDF8c2l0ZS53eTIpICsgKDF8cGl0X2lkKSwgZmFtaWx5ID0gYmlub21pYWwoKSwKICAgICAgICAgICAgY29udHJvbD1nbG1lckNvbnRyb2wob3B0aW1pemVyPSJib2J5cWEiLCBvcHRDdHJsPWxpc3QobWF4ZnVuPTEwMDAwMCkpLAogICAgICAgICAgICBkYXRhID0gc3Vic2V0KGRmLm0zKSk7IHN1bW1hcnkobTMpCmBgYAoKCmBgYHtyfQojIHZpZXcgdGVtcGVyYXR1cmUgc2xvcGVzIGZvciBlYWNoIHNleCB1c2luZyBlbXRyZW5kcwprYWJsZShtMy5lbXQpCmBgYAoKICAtICoqX1Jlc3VsdDpfKiogX0ZlbWFsZXMgYXJlIGxlc3MgYWN0aXZlIHRocm91Z2hvdXQgdGhlaXIgc3dhcm0gcGVyaW9kIGdlbmVyYWxseSwgYnV0IGZlbWFsZSBhY3Rpdml0eSBpbmNyZWFzZXMgd2l0aCB0ZW1wZXJhdHVyZS4gTWFsZSBhY3Rpdml0eSBpcyBsZXNzIGluZmx1ZW5jZWQgYnkgdGVtcGVyYXR1cmUuXwoKIyMgTW9kZWwgNCwgRmlndXJlIDNCOiBUaGVuLCB3aGF0IGFyZSB0aGUgZWZmZWN0cyBvZiB0aGUgbGVuZ3RoIG9mIHRoZSBhY3RpdmUgcGVyaW9kIGFuZCBzZXggb24gdGhlICoqX3Byb2JhYmlsaXR5IG9mIGxhc3QgZGV0ZWN0aW9uXyoqPwogIApgYGB7cn0KIyByZXNwb25zZTogMSA9IGxhc3QgZGF5IGFjdGl2ZSB8IDAgPSBkZXRlY3RlZCBvbiBzdWJzZXF1ZW50IG5pZ2h0cwptNCA9IGdsbWVyKGxhc3QuZGV0Mn5pc2Eqc2V4ICsgKDF8c2l0ZS53eTIpICsgKDF8cGl0X2lkKSwgZmFtaWx5PWJpbm9taWFsKCksCiAgICAgICAgICAgY29udHJvbD1nbG1lckNvbnRyb2wob3B0aW1pemVyPSJib2J5cWEiLCBvcHRDdHJsPWxpc3QobWF4ZnVuPTEwMDAwMCkpLAogICAgICAgICAgIGRhdGE9c3Vic2V0KGRmLm00KSk7c3VtbWFyeShtNCkKYGBgCgogIC0gKipfUmVzdWx0Ol8qKiBfUmVnYXJkbGVzcyBvZiBhcnJpdmFsIGRhdGUsIGZlbWFsZXMgc3BlbmQgZmV3ZXIgZGF5cyBzd2FybWluZyB0aGFuIG1hbGVzIGJlZm9yZSBlbmRpbmcgYWN0aXZpdHkuXwogIAojIyBNb2RlbCA1LCBGaWd1cmUgM0M6IExhc3QsIHdoYXQgYXJlIHRoZSBlZmZlY3RzIG9mIHRlbXBlcmF0dXJlIGFuZCBzZXggb24gdGhlICoqX3Byb2JhYmlsaXR5IG9mIGxhc3QgZGV0ZWN0aW9uP18qKgpgYGB7cn0KIyByZXNwb25zZTogMSA9IGxhc3QgZGF5IGFjdGl2ZSB8IDAgPSBkZXRlY3RlZCBvbiBzdWJzZXF1ZW50IG5pZ2h0cwptNSA9IGdsbWVyKGxhc3QuZGV0Mn5iYXQudGF2ZypzZXggKyAoMXxwaXRfaWQpICsgKDF8c2l0ZS53eTIpLCBmYW1pbHkgPSBiaW5vbWlhbCgpLAogICAgICAgICAgICBjb250cm9sPWdsbWVyQ29udHJvbChvcHRpbWl6ZXI9ImJvYnlxYSIsIG9wdEN0cmw9bGlzdChtYXhmdW49MTAwMDAwKSksCiAgICAgICAgICAgIGRhdGEgPSBzdWJzZXQoZGYubTUpKTsgc3VtbWFyeShtNSkKYGBgCi0gKipfUmVzdWx0Ol8qKiBfRmVtYWxlcyB3aWxsIGVuZCBhY3Rpdml0eSBvbiB3YXJtZXIgbmlnaHRzIGNvbXBhcmVkIHRvIG1hbGVzIHRoYXQgcmVtYWluIGFjdGl2ZSBldmVuIGFzIHRlbXBlcmF0dXJlcyBkZWNyZWFzZS5fCgojIyBQbG90IEZpZ3VyZSAzLiBDb250cmlidXRpb24gb2YgdGVtcGVyYXR1cmUgYW5kIGxlbmd0aCBvZiBzd2FybSBvbiBzZXgtc3BlY2lmaWMgYWN0aXZpdHkgKE1vZGVscyAzLTUpCmBgYHtyfQpwcmludChwLmZpZzMpCmBgYAo8ZGl2IGNsYXNzID0gImFsZXJ0IGFsZXJ0LWluZm8iIHJvbGU9ImFsZXJ0Ij4gIAogIC0gKipfS2V5IGZpbmRpbmdzXyoqOiBfRmVtYWxlIGFjdGl2aXR5IGlzIGdlbmVyYWxseSByZWR1Y2VkIGNvbXBhcmVkIHRvIG1hbGVzLiBGZW1hbGVzIGxpbWl0IG5pZ2h0bHkgYWN0aXZpdHkgdG8gb2NjdXIgYXQgd2FybWVyIHRlbXBlcmF0dXJlcywgYW5kIGZlbWFsZXMgYWxzbyBlbmQgYWN0aXZpdHkgZmV3ZXIgZGF5cyBhZnRlciBhcnJpdmluZyBhdCB0aGUgaGliZXJuYWN1bGEgYW5kIGF0IGhpZ2hlciB0ZW1wZXJhdHVyZXMuXwo8L2Rpdj4KCiMgSG93IGRvZXMgbWF0aW5nIHBoZW5vbG9neSBjb3JyZXNwb25kIHRvIHNlYXNvbmFsIGRpc2Vhc2UgZHluYW1pY3M/IChGaWd1cmUgNEEtQiwgTW9kZWxzIDYtOCkKCiMjIE1vZGVsIDY6IEZpcnN0LCBkbyAqKl9wYXRob2dlbiBsb2Fkc18qKiBkaWZmZXIgYmV0d2VlbiBzZXhlcyBkdXJpbmcgdGhlIGF1dHVtbiBtYXRpbmcgKGFjdGl2ZSkgc2Vhc29uPwpgYGB7cn0KbTYgPSBsbWVyKGxnZEwyfnNleCArICgxfHNpdGUud3kyKSwgCiAgICAgICAgICBkYXRhID0gZGYubTYpOyBzdW1tYXJ5KG02KSAKYGBgCiAgLSAqKl9SZXN1bHQ6XyoqIF9ObyBjbGVhciBpbmZlY3Rpb24gYmlhcyBkdXJpbmcgYXV0dW1uIG1hdGluZy5fCgojIyBNb2RlbCA3LCBGaWd1cmUgNEE6IFRoZW4sIGhvdyBkbyAqKl9wYXRob2dlbiBsb2Fkc18qKiBjaGFuZ2Ugc2Vhc29uYWxseSBieSBzZXg/CmBgYHtyfQptNyA9IGxtZXIobGdkTDJ+c2V4KnBkYXRlMiArICgxfHNpdGUud3kyKSwgCiAgICAgICAgICBkYXRhID0gZGYyKTsgc3VtbWFyeShtNykKYGBgCgogIC0gKipfUmVzdWx0XyoqOiBfRmVtYWxlcyBkZXZlbG9wIG1vcmUgc2V2ZXJlIGluZmVjdGlvbnMgYnkgZWFybHkgaGliZXJuYXRpb24uXwoKIyMgTW9kZWwgOCwgRmlndXJlIDRCOiBMYXN0LCBkb2VzIHRoZSBlbmQgZGF0ZSBvZiBhdXR1bW4gYWN0aXZpdHkgaW5mbHVlbmNlICoqX2luZmVjdGlvbnMgaW4gZWFybHkgaGliZXJuYXRpb25fKio/CmBgYHtyfQptOCA9IGxtZXIobGdkTDJ+bWVkLm1heGRhdGUgKyAoMXxzaXRlLnd5MiksIAogICAgICAgICAgIGNvbnRyb2w9bG1lckNvbnRyb2wob3B0aW1pemVyPSJib2J5cWEiLCBvcHRDdHJsPWxpc3QobWF4ZnVuPTEwMDAwMCkpLAogICAgICAgICAgIGRhdGEgPSBkZi5tOCk7IHN1bW1hcnkobTgpCmBgYAoKIyMgUGxvdCBGaWd1cmUgNC4gU2Vhc29uYWwgc2V4LWJpYXNlZCBpbmZlY3Rpb24gKEZpZ3VyZSA0QS1COyBNb2RlbHMgNy04KQoKYGBge3J9CnByaW50KHAuZmlnNCkKYGBgCjxkaXYgY2xhc3MgPSAiYWxlcnQgYWxlcnQtaW5mbyIgcm9sZT0iYWxlcnQiPiAgCiAgLSAqKl9LZXkgZmluZGluZ3M6XyoqIAogICAgLSBfRmVtYWxlLWJpYXNlZCBpbmZlY3Rpb24gb25seSBhcmlzZXMgZm9sbG93aW5nIHRoZSBhdXR1bW4gbWF0aW5nIHBlcmlvZC5fIAogICAgLSBfSW5mZWN0aW9ucyBhcmUgbW9yZSBzZXZlcmUgb24gYmF0cyBhc3NvY2lhdGVkIHdpdGggc2l0ZXMgd2hlcmUgYmF0cyBlbmRlZCBhY3Rpdml0eSBlYXJsaWVyLl8KPC9kaXY+CgogIAojIFN1cHBsZW1lbnRhbCBhbmFseXNlcyBhbmQgZmlndXJlcwoKIyMgU3VwcCBGaWd1cmUgMTogVmlzdWFsaXplIGNvbXBsZXRlIGRhdGFzZXQKCmBgYHtyfQpwcmludChwLnNmMSkKYGBgCiAgLSBHZW5lcmFsbHksIG1hbGVzIGFyZSBtb3JlIGFjdGl2ZSB0aGFuIGZlbWFsZXMgdGhyb3VnaG91dCBhdXR1bW4gbWF0aW5nLCBmdWxseSBlbmNvbXBhc3NpbmcgZmVtYWxlIHN3YXJtIGFjdGl2aXR5LgoKIyMgVXNlIGJhbGFuY2VkIGRhdGFzZXQgdG8gc3VwcG9ydCBzZXgtc3BlY2lmaWMgZWZmZWN0cwoKIyMjIFN1cHAgTW9kZWwgMSwgU3VwcCBGaWcgMkE6IFRlc3QgdGVtcGVyYXR1cmUtZGVwZW5kZW5jZSBvZiBiYXQgYWN0aXZpdHkKCmBgYHtyfQpzbTEgPSBnbG1lcihkZXRlY3QyfmJhdC50YXZnKnNleCArICgxfHNpdGUud3kyKSArICgxfHBpdF9pZCksIGZhbWlseSA9IGJpbm9taWFsKCksCiAgICAgICAgICAgICBjb250cm9sPWdsbWVyQ29udHJvbChvcHRpbWl6ZXI9ImJvYnlxYSIsIG9wdEN0cmw9bGlzdChtYXhmdW49MTAwMDAwKSksCiAgICAgICAgICAgICBkYXRhID0gc3Vic2V0KGRmMS50cnUpKTsgc3VtbWFyeShzbTEpCmBgYAoKYGBge3J9CiN2aWV3IHRlbXBlcmF0dXJlIHRyZW5kcyBieSBzZXgKa2FibGUoc20xLmVtdCkKYGBgCgojIyMgU3VwcCBNb2RlbCAyLCBTdXBwIEZpZyAyQjogRGlmZmVyZW5jZXMgaW4gbmlnaHRseSB0ZW1wZXJhdHVyZSBieSBzZXggYW5kIGFjdGl2aXR5CgpgYGB7cn0Kc20yID0gbG1lcihiYXQudGF2Z35kZXRlY3QyKnNleCArICgxfHBpdF9pZCkgKyAoMXxzaXRlLnd5MiksIAogICAgICAgICAgIGRhdGEgPSBzdWJzZXQoZGYxLnRydSkpOyBzdW1tYXJ5KHNtMikKYGBgCiAgCiMjIyBQbG90IFN1cHAgRmlndXJlIDIuIEVmZmVjdHMgb2YgdGVtcGVyYXR1cmUgd2l0aCBiYWxhbmNlZCBvYnNlcnZhdGlvbnMgYmV0d2VlbiBzZXhlcwpgYGB7cn0KcHJpbnQocC5zZjIpCmBgYAogIC0gKipfUmVzdWx0cyBjb25zaXN0ZW50IHdpdGggTW9kZWwgMyBpbiBGaWd1cmUgM0EsIHVzaW5nIHRoZSBiYWxhbmNlZCBkYXRhc2V0IGJldHdlZW4gc2V4ZXM6XyoqIAogICAgLSBfRmVtYWxlcyBpbmNyZWFzZSBuaWdodGx5IGFjdGl2aXR5IHdpdGggdGVtcGVyYXR1cmUgY29tcGFyZWQgdG8gbWFsZXMgd2hvc2UgYWN0aXZpdHkgaXMgcmVsYXRpdmVseSB1bmFmZmVjdGVkIGJ5IHRlbXBlcmF0dXJlLl8gCiAgICAtIF9UaGUgZGlmZmVyZW5jZSBiZXR3ZWVuIG1lYW4gdGVtcGVyYXR1cmUgb24gbmlnaHRzIHdoZW4gYmF0cyB3ZXJlIGFjdGl2ZSBjb21wYXJlZCB0byBuaWdodCB3aGVuIGJhdHMgd2VyZSB1bmRldGVjdGVkIHdhcyBncmVhdGVyIGZvciBmZW1hbGVzIHRoYW4gbWFsZXMuIEZlbWFsZXMgY29uY2VudHJhdGUgdGhlaXIgYWN0aXZpdHkgdG8gd2FybWVyIG5pZ2h0cyBjb21wYXJlZCB0byBtYWxlcyB0aGF0IGRvIG5vdC5fCgojIyBEb2VzIGJvZHkgY29uZGl0aW9uIGFmZmVjdCBuaWdodGx5IGFjdGl2aXR5IChTdXBwIEZpZ3VyZSAzOyBTdXBwIE1vZGVsIDMpPwoKIyMjIFRlc3Qgd2hldGhlciB0ZW1wZXJhdHVyZSBhbmQgbWFzcyBkaWZmZXJlbnRseSBpbmZsdWVuY2UgYWN0aXZpdHkgYmV0d2VlbiBzZXhlcwoKYGBge3J9CnNtMyA8LSBnbG1tVE1CKGRldGVjdDJ+c2V4Km1hc3MqdGVtcCArICgxfHBpdF9pZCksCiAgICAgICAgICAgICAgIGZhbWlseSA9IGJpbm9taWFsKCksCiAgICAgICAgICAgICAgIGRhdGEgPSBzdWJzZXQoZGYuc20zKSk7IHN1bW1hcnkoc20zKQpgYGAKCiMjIyBQbG90IFN1cHAgRmlndXJlIDM6IEVmZmVjdHMgb2YgbWFzcyBhbmQgdGVtcGVyYXR1cmUgb24gbmlnaHRseSBhY3Rpdml0eQoKYGBge3J9CnByaW50KHAuc2YzKQpgYGAKCiAgLSAqKl9SZXN1bHRzOl8qKiAKICAgICAgLSBfRmVtYWxlcyB3aXRoIGxvd2VyIGJvZHkgY29uZGl0aW9uIChpZSBib2R5IG1hc3MpIGNvbmNlbnRyYXRlZCBhY3Rpdml0eSB0byB3YXJtIG5pZ2h0cywgY29tcGFyZWQgdG8gZmVtYWxlcyB3aXRoIGhpZ2hlciBib2R5IGNvbmRpdGlvbi5fCiAgICAgIC0gX1RoZSB0ZW1wZXJhdHVyZXMgYXQgd2hpY2ggbWFsZXMgd2VyZSBhY3RpdmUgd2FzIG5vdCBhZmZlY3RlZCBieSB0aGVpciBib2R5IGNvbmRpdGlvbi5fCiAgICAgIC0gX0ZlbWFsZXMgYWRqdXN0IG5pZ2h0bHkgYWN0aXZpdHkgYWNjb3JkaW5nIHRvIHRoZWlyIGJvZHkgY29uZGl0aW9uIGJ1dCBtYWxlcyBkbyBub3QsIHN1Z2dlc3RpbmcgdGhlIHNleGVzIGFyZSBsaWtlbHkgYnVkZ2V0aW5nIHRoZWlyIGVuZXJneSBkaWZmZXJlbnRseS5fCgoKIyMgQWRkaXRpb25hbCBhbmFseXNpcyBzaG93aW5nIGFzc29jaWF0aW9uIGJldHdlZW4gYmV0d2VlbiBtYXRpbmcgcGhlbm9sb2d5IGFuZCBkaXNlYXNlIChTdXBwIEZpZ3VyZSA0OyBTdXBwIEZpZ3VyZSA0KQoKIyMjIFRlc3Qgd2hldGhlciBtZWRpYW4gZW5kIGRhdGUgb2YgYXV0dW1uIGFjdGl2aXR5IGF0IGEgc2l0ZSBpbmZsdWVuY2VkIG1lYW4gcGF0aG9nZW4gbG9hZHMgaW4gZWFybHkgaGliZXJuYXRpb24KCmBgYHtyfQpzbTQgPSBsbWVyKG1lYW4ubGdkTH5tZWQubWF4ZGF0ZSArICgxfHNpdGUud3kyKSwgZGF0YSA9IHN1YnNldChkZjIuc3VtKSk7IHN1bW1hcnkoc200KQpgYGAKCgojIyMgUGxvdCBTdXBwIEZpZ3VyZSA0LiBFZmZlY3Qgb2YgbWVkaWFuIGVuZCBkYXRlIG9mIGF1dHVtbiBhY3Rpdml0eSBvbiBtZWFuIGVhcmx5IGhpYmVybmF0aW9uIGluZmVjdGlvbnMgYXQgc2l0ZS4KCmBgYHtyfQpwcmludChwLnNmNCkKYGBgCiAgLSAqKl9SZXN1bHRzOiBDb25zaXN0ZW50IHdpdGggRmlndXJlIDRCLCBpbmZlY3Rpb24gc2V2ZXJpdHkgaW4gZWFybHkgaGliZXJuYXRpb24gaXMgbW9yZSBzZXZlcmUgYXQgc2l0ZXMgd2hlcmUgYXV0dW1uIGFjdGl2aXR5IGVuZGVkIHNvb25lci5fKiogCgojIyBIb3cgbWFueSBuaWdodHMgb2YgYWN0aXZpdHkgZG8gaW5kaXZpZHVhbHMgYWNjdW11bGF0ZSB0aHJvdWdob3V0IHRoZSBsZW5ndGggb2YgdGhlaXIgYWN0aXZlIHBlcmlvZD8gCgojIyMgUGxvdCBTdXBwIEZpZ3VyZSA1OiBWaXN1YWxpemUgY3VtdWxhdGl2ZSBudW1iZXIgb2YgYWN0aXZlIG5pZ2h0cyBieSBzZXgKCmBgYHtyfQpwcmludChwLnNmNSkKYGBgCgojIyMgQ2FsY3VsYXRlIGF2ZXJhZ2UgbnVtYmVyIG9mIG5pZ2h0cyBkZXRlY3RlZCBhbW9uZyBpbmRpdmlkdWFscyBvZiBlYWNoIHNleAoKYGBge3J9CmthYmxlKG4uZGF0ZXMuc3VtKQpgYGAKCiAgLSBNYWxlcyBjb250aW51ZSB0byBiZSBhY3RpdmUgd2l0aG91dCBhbnkgY2xlYXIgdGhyZXNob2xkIGluIHRoZSBudW1iZXIgb2YgZGF0ZXMgYWN0aXZlLCB3aGVyZWFzIGZlbWFsZXMgY2Vhc2UgYWN0aXZpdHkgYWZ0ZXIgYW4gYXZlcmFnZSBvZiA0IG5pZ2h0cy4KICAKIyMgTW9kZWwgY29tcGFyaXNvbnMgYW5kIHZhbGlkYXRpb24KCiMjIyBDb21wYXJlIGFjdGl2aXR5IHBlcmlvZCBtb2RlbHMgYXMgYWRkaXRpdmUgb3IgaW50ZXJhY3RpdmUgd2l0aCBzaXRlLXllYXIgYW5kIHNleCBhcyBmaXhlZCBlZmZlY3RzIChTdXBwIFRhYmxlIDIpCgogRmlyc3QsIGJlZ2luIGRhdGUgb2YgYWN0aXZpdHkgbW9kZWxzCmBgYHtyfQptMWIuYWRkID0gbG1lcihtaW5kYXRlMn5zaXRlLnd5MiArIHNleCArICgxfHBpdF9pZCksIGRhdGEgPSBkZi5tMSkKbTFiLmludCA9IGxtZXIobWluZGF0ZTJ+c2l0ZS53eTIgKiBzZXggKyAoMXxwaXRfaWQpLCBkYXRhID0gZGYubTEpCmBgYAoKCiAtIENvbXBhcmUgdXNpbmcgQUlDCmBgYHtyfQprYWJsZShBSUMobTFiLmFkZCxtMWIuaW50KSkKYGBgCgoKIFRoZW4sIGVuZCBkYXRlIG9mIGFjdGl2aXR5IG1vZGVscwpgYGB7cn0KbTJiLmFkZCA9IGxtZXIobWF4ZGF0ZTJ+c2l0ZS53eTIgKyBzZXggKyAoMXxwaXRfaWQpLCBkYXRhID0gZGYubTIpCm0yYi5pbnQgPSBsbWVyKG1heGRhdGUyfnNpdGUud3kyICogc2V4ICsgKDF8cGl0X2lkKSwgZGF0YSA9IGRmLm0yKQpgYGAKCiAgLSBDb21wcmUgdXNpbmcgQUlDCmBgYHtyfQprYWJsZShBSUMobTJiLmFkZCxtMmIuaW50KSkKYGBgCgo8ZGl2IGNsYXNzID0gImFsZXJ0IGFsZXJ0LWluZm8iIHJvbGU9ImFsZXJ0Ij4gIAogIC0gKipfS2V5IHRha2Vhd2F5Ol8qKgogICAgLSBfU2V4IGlzIG1vcmUgc3VwcG9ydGVkIGFzIGFuIGFkZGl0aXZlIHJhdGhlciB0aGFuIGludGVyYWN0aXZlIGVmZmVjdCBpbiBtb2RlbHMgaW5jbHVkaW5nIHNpdGUteWVhciBhcyBmaXhlZCBlZmZlY3RzLl8KICAgIC0gX1RoZSBkYXRlcyB0aGF0IGFjdGl2aXR5IHN0YXJ0ZWQgYW5kIGVuZGVkIHZhcmllZCBhbW9uZyBzaXRlcyBhbmQgeWVhcnMsIGJ1dCBmZW1hbGVzIGNvbnNpc3RlbnRseSB3ZXJlIGxhdGVyIHRvIHN0YXJ0IGFuZCBlYXJsaWVyIHRvIGVuZC5fCjwvZGl2PgoKCiMjIyBWaWV3IHJlc3VsdHMgb2YgbnVsbCBtb2RlbCBjb21wYXJpc29ucyBhbmQgQVVDIGVzdGltYXRlcyBmcm9tIGstZm9sZCBjcm9zcyB2YWxpZGF0aW9uIG9mIGJpbm9taWFsIG1vZGVscwoKYGBge3J9CmthYmxlKG51bGxjb21wc19mdWxsKQpgYGAKCiAgLSBBbGwgQUlDcyBvZiByZXBvcnRlZCBtb2RlbHMgYXJlID4yIHNjb3JlcyBiZWxvdyB0aGUgbnVsbCBtb2RlbHMgKHN0cnVjdHVyZWQgd2l0aCByZXNwb25zZX4xKSwgaW5kaWNhdGluZyBpbXByb3ZlbWVudCBvdmVyIHRoZSBudWxsLiAKCiMjIyBWaWV3IEFVQyBzY29yZXMgZGVyaXZlZCBmcm9tIGstZm9sZCBjcm9zcyB2YWxpZGF0aW9uIG9mIHRoZSBsb2dpc3RpYyBtb2RlbHMgcmVwb3J0ZWQgaW4gdGhlIG1haW4gcmVzdWx0cy4KYGBge3J9CmthYmxlKGtmb2xkX2Z1bGwpCmBgYAogICAtIFJlc3VsdHMgaW5kaWNhdGUgdGhhdCA+NzAlLCA3NSUsIGFuZCA4MyUgb2Ygb3VyIHRlc3QgZGF0YSB3YXMgc3VjY2Vzc2Z1bGx5IHByZWRpY3RlZCBieSB0aGUgdHJhaW5pbmcgbW9kZWxzIG9mIE1vZGVsIDMsIE1vZGVsIDQsIGFuZCBNb2RlbCA1LCByZXNwZWN0aXZlbHkuCgojIyBWaWV3IHNhbXBsaW5nIGRpc3RyaWJ1dGlvbnMgYnkgc2l0ZSwgeWVhciBhbmQgc2V4IHdoZXJlIGFwcGxpY2FibGUgYnkgbW9kZWwKYGBge3IsIGVjaG89RkFMU0UscmVzdWx0cz0naGlkZSd9CnByaW50KHNhbXBsZXNpemVzX2Z1bGwpCmBgYAoK